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Antiferromagnetic Heisenberg Spin Chain of a Few Cold Atoms in a One-Dimensional Trap 
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We report on the deterministic preparation of antiferromagnetic Heisenberg spin chains consisting of up to 
four fermionic atoms in a one-dimensional trap. These chains are stabilized by strong repulsive interactions 
between the two spin components without the need for an external periodic potential. We independently char¬ 
acterize the spin configuration of the chains by measuring the spin orientation of the outermost particle in the 
trap and by projecting the spatial wave function of one spin component on single-particle trap levels. Our 
results are in good agreement with a spin-chain model for fermionized particles and with numerically exact 
diagonalizations of the full few-fermion system. 


The high control and tunability of ultracold atomic systems 
offer the fascinating possibility to simulate quantum mag¬ 
netism [1], a topic of fundamental importance in condensed 
matter physics [2], Systems of spin-1/2 fermions with antifer¬ 
romagnetic (AFM) correlations are of particular interest due 
to the observation of high-temperature superconductivity in 
cuprates with AFM correlations [3], The experimental imple¬ 
mentation of the necessary exchange couplings is usually real¬ 
ized by superexchange processes of neighboring atoms in the 
Mott-insulating state of a deep optical lattice. Superexchange 
couplings were measured in both bosonic [4] and fermionic 
double-well systems [5] and short-range AFM correlations 
of fermionic atoms were detected in various lattice geome¬ 
tries [6-8]. Furthermore, superexchange processes were used 
to study the dynamics of spin impurities above the ferromag¬ 
netic (FM) ground state of bosons in the Mott-insulating state 
of a one-dimensional (ID) lattice [9], Bosonic atoms were 
also used to simulate AFM Ising spin chains in a tilted optical 
lattice [10, 11]. However, the AFM ground state of spin-1/2 
fermions in a deep optical lattice has so far not been realized 
due to the very low energy scale associated with the superex¬ 
change coupling. 

This problem can be circumvented in ID systems, where 
quantum magnetism can be simulated without an optical lat¬ 
tice [12-14]. In the regime of strong interactions, the spa¬ 
tial wave function of both fermions [15] and bosons [16- 
18] can be mapped on the wave function of spinless nonin¬ 
teracting fermions [Fig. 1(a)]. In this so called fermioniza- 
tion limit, the strong interactions lead to the formation of a 
Wigner-crystal-like state [19-21], which has a highly degen¬ 
erate ground state when the particles have multiple internal 
degrees of freedom [Fig. 1(b)] [20-23], Close to the limit of 
fermionization, the structure of the quasi-degenerate ground- 
state multiplet [24-33] is determined by an effective Suther¬ 
land spin-chain Hamiltonian, which for two-component sys¬ 
tems becomes a Heisenberg model [12, 19, 21, 29, 32-34], 

In this paper, we report on the realization of Heisenberg 
spin chains of A/- spin-up and N\, spin-down particles with 
(Nf, iVj,) = (2,1), (3,1), and (2, 2). We show that under an 
adiabatic change of the interaction strength the noninteracting 
ground states of these systems evolve into the respective AFM 
states in the limit of infinitely strong repulsion [24, 28], We 


identify the AFM states by two independent measurements. 
First, we use a tunneling technique to measure the spin orien¬ 
tation of the outermost particle of the spin chain. Second, we 
probe the spatial wave function of the spin-down atom in the 
(2,1) and (3,1) system by projecting it on single-particle trap 
levels. 

In our experiments, we realize a spin-1/2 system 
by trapping ultracold b Li atoms in an elongated op¬ 
tical dipole trap [35, 39] in their two lowest hyper- 
fine states |t) = | j = 1/2,mj = —1/2; I = 1 ,m/ = 0) and 
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FIG. 1: Heisenberg spin chain of three fermions, (a) Sketch of two 
spin-up and one spin-down atom with diverging ID coupling con¬ 
stant (pm = ±oo) in a harmonic trap. If the relative spatial wave 
function of two distinguishable fermions is symmetric, the strong in¬ 
teractions induce a cusp in the relative wave function of the two par¬ 
ticles (left side). This causes them to separate like identical fermions 
(right side). In this fermionization limit the system forms a Wigner- 
crystal-like state with fixed ordering of the particles, (b) Single¬ 
particle contributions to the total (gray), the spin-up (green), and the 
spin-down density (blue) of two spin-up and one spin-down atom in 
the fermionization regime in a harmonic trap. Like in a Wigner crys¬ 
tal, the total densities of the ferromagnetic (left), the intermediate 
(middle), and the antiferromagnetic state (right) are identical, while 
their spin densities differ and are determined by a Heisenberg spin- 
chain Hamiltonian, (c) Densities of three particles before (left) and 
after (right) the tunneling of one atom with energy E out of a tilted 
trap. At fermionization, only the rightmost particle can leave the trap 
in the tunneling process. 
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I4-) = |j = 1 / 2 , = — 1 / 2 ;/ = 1 , rnj = 1 ). As the energy 

of the atoms is much smaller than the lowest transverse ex¬ 
citation energy in the trap, their dynamics are restricted to 
the longitudinal axis of the trap. In such a quasi-ID system, 
the interaction strength between ultracold atoms of opposite 
spin is determined by the ID coupling constant, gin, which 
diverges at a confinement-induced resonance (CIR) when the 
3D scattering length, 030 , approaches the harmonic oscillator 
length of the radial confinement [35, 40], We use a magnetic 
Feshbach resonance to control 0,30 and therefore are able to 
smoothly tune g id across the CIR. At the same time, scatter¬ 
ing between fermionic atoms of the same spin component is 
forbidden. Throughout this paper, g id will be given in units of 
a\\hhj\\, where a|| = yJhfmLj |j and W|| are the harmonic os¬ 
cillator length and the trap frequency in longitudinal direction 
and m is the mass of a 6 Li atom. 

We start our experiments by preparing a (2,1), (3,1), or 
( 2 , 2 ) system in the noninteracting many-particle ground state 
of the trap [35, 39]. By changing the magnetic offset field with 
a constant rate, we ramp the system into the fermionization 
regime close to the CIR (Fig. 2), where it forms a spin chain. 
Below the CIR, <?id is large and positive and the system is in 
the Tonks regime of strong repulsion [17, 18]. When crossing 
the CIR, g\n changes sign from +00 to —00 while the system 
continuously follows the so called upper branch [24] into the 
super-Tonks regime of strong attraction [15, 41, 42] (Fig. 2). 
In the super-Tonks regime, the system is in an excited state, 
which is metastable against decay into bound states. 

In a first set of measurements, we identify the states of the 
spin chains by probing the spin distributions in the trap. Here, 
we make use of the fact that in the fermionization regime the 
atoms become impenetrable and therefore their ordering along 
the longitudinal axis of the trap is fixed. This allows us to de¬ 
termine the spin orientation of the outermost particle in the 
trap in a tunneling measurement. To do this, we tilt the trap 
as shown in Fig. 1(c) and thereby allow atoms to tunnel out 
of the trap. We carefully adjust the trap parameters during the 
tunneling process, to let exactly one atom (for the ( 2 , 1 ) and 
the (3,1) systems) or two atoms (for the (2, 2) system) tun¬ 
nel [35]. Finally, we measure the number of spin-up atoms in 
the final state to determine the spin of the atoms that left the 
trap during the tunneling process [35]. We define spin-down 
tunneling as the process in which all spin-down atoms tunnel 
out. By repeating this measurement at different magnetic off¬ 
set fields, we deduce the probability of spin-down tunneling, 
P^(—1 /(?id), as a function of the inverse ID coupling constant 
as shown in Fig. 3. 

As shown in Fig. 1(b) for a (2,1) system in a harmonic trap, 
the different states of the spin chain can be uniquely identi¬ 
fied by their spin densities [ 12 ], specifically by the probabil¬ 
ity of the outermost spin to point downwards. Since in the 
fermionization regime the ordering of the atoms in the trap is 
fixed only the outermost atom can escape during the tunnel¬ 
ing process. Exactly at the CIR, the probability of spin-down 
tunneling should therefore directly reveal the state of the spin 
chain [12, 29, 33], Away from resonance, the probability of 



FIG. 2: Energies in spin-chain regime. Eigenenergies of two (green) 
and three (red) strongly interacting spin-1/2 fermions in a ID har¬ 
monic trap as a function of the interaction strength. In the Tonks 
regime, the antiferromagnetic states are the ground states of each 
multiplet, while the ferromagnetic states have the highest energies. In 
the super-Tonks regime, the ordering of the energy levels is inverted. 
Close to the confinement-induced resonance (CIR), the energy shifts 
are linear in — 1/gm and can be determined by a Heisenberg spin- 
chain Hamiltonian. The system is initially prepared in the noninter¬ 
acting ground state of the three-particle system at — 1/gm = — 00 , 
which evolves, for increasing — 1/<?id, into the antiferromagnetic 
state around the CIR (red solid line). During a ramp across the CIR, 
the system stays in the antiferromagnetic state, since all eigenstates 
of the system are decoupled. The blue arrows indicate the predom¬ 
inant channels for the tunneling of one atom below (left) and above 
(right) the fermionization regime. 


spin-down tunneling is also influenced by the energy of the 
final in-trap states, favoring final states with lower energy as 
indicated by the blue arrows in Fig. 2. To identify the spin 
states throughout the entire spin-chain regime, we compare 
our data to the results of a tunneling model, which in the fol¬ 
lowing section is explained for a ( 2 , 1 ) system. 

In our tunneling model, the initial states are eigenstates of a 
Heisenberg spin-chain Hamiltonian [35], where the exchange 
couplings Ji between neighboring spins depend on the trap 
geometry and on the inverse ID coupling constant [12], For 
the ( 2 , 1 ) system with repulsive interactions and a symmetric 
trap (Ji = J 2 > 0 ), these eigenstates are the antiferromag¬ 
netic (AFM) ground state, the intermediate (IM) state, and 
the ferromagnetic (FM) state as shown in Fig. 1(b). During 
the tunneling process the trap is tilted as shown in Fig. 1(c) 
and therefore the density is not symmetric. Hence, the ex¬ 
change couplings are not identical anymore (| Ji I > \J2\), 
which leads to a coherent mixing of the AFM and IM state 
during the tunneling process [35]. We calculate a probability 
of approximately 8 % for the rightmost spin of the AFM state 
in the tilted trap to point downwards. This is in good agree¬ 
ment with the blue data points in Fig. 3(a) that cross the CIR 
at P i « 10%. 

Away from the CIR, the eigenstates of both the three- 
particle and the two-particle spin chains are nondegenerate 
(Fig. 2). In this case, the energies of the initial three-particle 
state | i) and the final two-particle state |/) involved in the 
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FIG. 3: Probing the spin distribution. Tunneling probabilities of the 
spin-down atom in a (2,1) system (a) and a (3,1) system (b) and tun¬ 
neling probability of both spin-down atoms in a (2, 2) system (c) as 
a function of the interaction strength. The red lines are the solutions 
of a tunneling model for the antiferromagnetic (solid), the ferromag¬ 
netic (dashed), and the intermediate state (dotted). The gray points 
in (a) indicate a narrow resonance between the antiferromagnetic and 
the intermediate state of the (2,1) system close to — 1/giD = 0. Er¬ 
ror bars denote the lcr statistical uncertainties. 


tunneling process are important, since their difference deter¬ 
mines the energy E of the tunneling particle. The tunneling 
rate of the particle that leaves the trap is strongly affected by 
its energy and can be calculated as 

T iJ( x\{i\f,t)\ 2 Ee~ 2 ^ E \ (1) 


where | f,t) = \f) < 8 > \ t) with 1 t) indicating the spin orientation 
of the tunneling particle. The tunneling parameter 7 is deter¬ 
mined by means of a WKB calculation [35]. The probability 
to tunnel from state |i) to state \f) is given by 


Pi,f = 


T, 


ij 


(E f .T ilf .y 


( 2 ) 


where the sum is over all possible final states | f). 

Using Eq. (2), we calculate the probabilities of tun¬ 

neling into the spin-polarized final state [red lines in Fig. 3(a)], 
which is equivalent to the probability of spin-down tunnel¬ 
ing (P^). Far below the CIR, the energy dependent term 
Ee ~ 27 dominates the outcome of the tunneling rates (Eq. 1). 
Therefore, tunneling into the AFM two-particle ground state 
(|t, 4-) — |t, t))/v^ i s strongly favored if its spin overlap to 
the initial state is not zero. This leads to a limiting value of 
Pi = 0 for initial AFM and IM states. Above the resonance, 
the energy ordering of the two-particle FM and AFM states 
is reversed and tunneling into the FM states is predominant 
(Fig. 2) [44]. Here, 7/ is determined by the ratio of the spin 
overlaps between the first two spins of the initial states and the 
FM two-particle states |f,t) and + |!, t))/\/2- 


The comparison of the theoretically predicted Pi of the 
AFM state in the tilted trap [red solid line in Fig 3(a)] with the 
experimental data (blue points) shows good agreement, while 
the FM (red dashed line) and IM (red dotted line) states are 
clearly excluded. We therefore conclude that before tunneling 
both below and above the CIR the system is in the AFM state. 
The gray points at — 1/<?id ~ 0 indicate a narrow resonance 
effect that couples the AFM state to the IM state of the spin 
chain. Since this resonance is accompanied by strongly en¬ 
hanced three-body losses [35], we suspect it to be caused by a 
coupling of the AFM and the IM states via a molecular state 
with center-of-mass excitation. The coupling to such molec¬ 
ular states is strongly enhanced by the anharmonicity of our 
tilted trap [43]. 

For the AFM state of the (3,1) system, a similar calcula¬ 
tion predicts // « 1% on resonance and a saturation value of 
Pi ss 75% deep in the super-Tonks regime [35]. As shown 
in Fig 3(b), the general trend of our measurements agrees 
with this prediction for the AFM state, but in the super-Tonks 
regime, there is a significant deviation. The reason for this de¬ 
viation is that the calculation assumes an adiabatic lowering 
of the potential barrier. As a result, the tunneling energies of 
all tunneling channels are always well below the barrier max¬ 
imum. We believe that this condition is not fulfilled for the 
(3,1) system in the super-Tonks regime, where an especially 
low potential barrier was used for the tunneling measurement. 
Indeed, if we model a nonadiabatic lowering of the potential 
barrier, the contribution from tunneling into the IM state re¬ 
duces Pi to values that are compatible with the experimental 
results [35], In order to study the spin configuration of the 
balanced (2, 2) system, we adapt the previous procedure and 
let two atoms tunnel out of the trap. Here, Pi is defined as 
the probability to end up in state |t, t), where both spin-down 
atoms tunneled out of the trap. Again, the predicted Pi « 4% 
on resonance and the limiting value of Pj_ « 33.3% in the 
super-Tonks regime are in good agreement with the experi¬ 
ment as shown in Fig. 3(c). 

To independently confirm the results of our measurement 
of the spin distribution, we perform a second set of measure¬ 
ments that directly probes the spatial wave function of the sys¬ 
tem. As shown in Fig. 1(a), the relative spatial wave function 
between identical spins always exhibits a smooth zero cross¬ 
ing, while between distinguishable spins with strong interac¬ 
tions it can exhibit a cusp. The cusps lead to occupancies of 
high-energy trap levels, while the zero crossings require only 
the occupation of the lowest trap levels. In general, the more 
symmetric the spatial wave function of a state is, the more 
cusps it will contain. Therefore, the occupation-number dis¬ 
tribution on single-particle trap levels directly reveals the spin 
configuration of the system. 

To probe this distribution, we prepare an interacting (2,1) 
or (3,1) system and remove all atoms of the spin-up com¬ 
ponent from the trap with a short pulse of light. The light 
is er~ polarized and resonant to the D2 transition of the 
spin-up atoms (|f) = \j = 1/2, nrij = —1/2; / = 1, to j = 0) 
to |j = 3/2,mj = —3/2; I = 1 ,mj = 0)). We confirm that 
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FIG. 4: Probing the spatial wave function. Occupation-number dis¬ 
tribution of the spin-down atom on single-particle trap levels for an 
initial (a) (2,1) and (b) (3,1) system. The red and gray symbols 
show theoretical predictions for the antiferromagnetic state (red cir¬ 
cles), the ferromagnetic state (gray squares), and intermediate states 
(gray diamonds). Both measurements (blue points) were made in the 
super-Tonks regime [— 1/gm = 0.586 ± 0.014 for the (2,1) system 
and — 1/t/io = 0.536 ± 0.013 for the (3,1) system] to show that 
both systems stay in the respective antiferromagnetic states through¬ 
out the regime of fermionization. Error bars denote the la statistical 
uncertainty. 


within our experimental fidelity all spin-up atoms are removed 
from the trap by the light pulse, while only 3 % of the pop¬ 
ulation of spin-down atoms is lost. With 15/iS the duration 
of the light pulse is significantly shorter than the inverse lon¬ 
gitudinal trap frequency of approximately 100//s, which sets 
the timescale of redistribution along the spin chain. This pro¬ 
cess therefore projects the spin-down component of the wave 
function of the interacting (Nf, l)-particle system on single¬ 
particle trap levels. Finally, we measure the mean occupancies 
on the single-particle trap levels [35]. In Fig. 4 we compare 
the mean occupancies of the spin-down atom for the (2,1) and 
the (3,1) systems in the super-Tonks regime with the theoret¬ 
ical prediction that we obtained by numerically diagonalizing 
the many-body Hamiltonian for these systems. The compar¬ 
ison shows that both systems are in the AFM spin state and 
thereby confirms that our systems follow this state throughout 
the fermionization regime. 

In conclusion, we have prepared antiferromagnetic Heisen¬ 
berg spin chains of up to four atoms in a one-dimensional trap 
and independently probed the spin distributions and spatial 
wave functions of the systems. This constitutes a direct 
observation of quantum magnetism beyond two-particle 
correlations in a system of ultracold fermionic atoms. By 
using the methods developed in Ref. [5], multiple spin chains 
can be realized and coupled, which offers a new approach to 
studying two and three-dimensional quantum magnetism. 
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Supplemental material 

I. Trapping potential 

The optical dipole trap (ODT) is created by the focus of a single far red-detuned Gaussian laser beam. For a given laser power 
P, the potential along the propagation direction of the light, which we denote the longitudinal axis, can be written as: 

KpticalM = PV 0 (l - 1 + {z/zr)2 )- (SI) 

Here, Vo = 3.326 /iK is the initial depth of the optical potential at a laser power of Pq = (265 ± 27) p.W, p = P/P 0 is the 

trap depth parameter and zr = is the Rayleigh range of the optical trapping beam with minimal waist wq = 1.838 /jm and 
wavelength A = 1064 nm. By changing the laser power P, we can adjust the trap depth during our experiments. In a harmonic 
approximation the trap frequencies of the dipole trap are W|| = 2n v /p (1.488 ± 0.014) kHz along the longitudinal direction and 
lo± = 2n^/p (14.22 ± 0.35) kHz in the radial direction, resulting in an aspect ratio of rj ss 10. 

To allow atoms to tunnel out of the trap, we create a potential barrier by superimposing the optical dipole trap with a linear 
magnetic gradient along the longitudinal axis [See Fig. 1(c) of main text]. This adds the magnetic part 

(Magnetic (^) = krn ^ Z (S2) 

to the potential, where //„, is the magnetic moment of the atoms and B' = 18.92 G/cm is the strength of the magnetic field 
gradient. For detailed information on the determination of the potential parameters, see Ref. [36], 


II. Spilling of atoms from the trap 

In several parts of our experimental sequence, we tilt the potential of the optical dipole trap (ODT) to spill atoms from the trap 
[Fig. 1(c) of the main text]. To do this, we apply a magnetic field gradient of B' = 18.92 G/cm, which adds a linear potential 
along the longitudinal axis of the ODT (Section I). The strength of this linear magnetic potential depends on the magnetic 
moment of the atoms in the trap, which is shown in Figure SI (a) as a function of the magnetic offset field. By tilting the ODT, 
we create a potential barrier that separates the in-trap states from the continuum of states outside of the trap. We reduce the 
height of the potential barrier by lowering the optical power in the ODT (p < 0). Atoms now leave the trap with rates that 
strongly depend on their energy. This allows us to remove atoms with energies above or marginally below the potential barrier 
from the trap while keeping atoms with significantly lower energies. After a certain time we ramp the optical power in the ODT 
back up to its initial value and thereby stop the spilling process. 

When performing this spilling technique in the Paschen-Back regime, the magnetic moments of spin-up and spin-down atoms 
are approximately equal [Fig. Sl(b)]. Therefore, the potential has the same shape for the spin-up and spin-down atoms. All 
tunneling measurements of the main text are performed deep in the Paschen-Back regime (> 725 G), where the magnetic 
moments of spin-up and spin-down atoms differ by less than 0.15 %. We use the spilling technique to let exactly one atom (for 
the (2,1) and the (3,1) system) or two atoms (for the (2,2) system) tunnel out of the trap (Section VII). 


III. Preparation of initial states and measurement of atom numbers 
Initial state preparation 

To prepare the initial few-particle systems, we start with a degenerate gas of about 600 noninteracting 6 Li atoms in their two 
lowest hyperfine states in the ODT. At a magnetic offset field of 527 G, where the atomic sample is noninteracting and the 
magnetic moments of spin-up and spin-down atoms are approximately equal, we spill all atoms above a certain trap level 
from the trap (Section II). This allows us to prepare spin-balanced ground-state systems of up to 10 atoms with a fidelity of 
~ 90 % [39], 

To prepare spin-imbalanced systems, we first prepare a spin-balanced system and then perform a second spilling process at 
a magnetic offset field of 40 G. At this field, the magnetic moment of the spin-up atoms (F = 1/2, m f = —1/2) is negligible 
[See inset of figure S 1(a)], which allows us to selectively spill spin-down atoms (F = 1/2, rriF = +1/2) from the trap [39], 
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FIG. SI: Magnetic moment of spin-up and spin-down atoms (a) and difference of the magnetic moments of spin-up and spin-down atoms as a 
function of the magnetic field. 


Atom number measurement 

In a single experimental run, we can either measure the total atom number N or the number of spin-up atoms TVf in the ODT. To 
measure N, we release the atoms from the ODT, recapture them in a magneto-optical trap and measure their fluorescence signal. 
To measure N^, we selectively spill all spin-down atoms from the trap and then measure the number of remaining atoms. 

A more detailed description of the preparation of both balanced and imbalanced few-fermion systems and of the detection of 
atom numbers in our experiment can be found in Ref. [39]. 


IV. Calculation of the ID coupling constant 

In a gas of ultracold 6 Li atoms, the interactions between the atoms are well described by contact s-wave scattering. In a one¬ 
dimensional system, the interaction potential between two atoms in different hyperfine states can be written as: 


VLu(21 - z 2 ) = gm S(zi - z 2 ), 


(S3) 


where Z\ and z 2 are the positions of the two atoms, g lD is the ID coupling constant, and S is the delta function [40], Contact 
interactions between identical fermions are forbidden. 

Starting from the 3D s-wave scattering length <230 [37], we calculate g id as: 


2fi, 2 d3D 1 

ma\ 1 - Ca 3 v/\/2a±_ ’ 


(S4) 


where m is the mass of a 6 Li atom, aj_ = y is the harmonic oscillator length of the radial confinement and 

C = —C(|) ~ 1-46 with the Riemann zeta function £ [40]. Figure S2 shows the calculated ID coupling constant <j\n between 
two atoms in the two lowest hyperfine states of 6 Li, which are trapped in our optical dipole trap with trap depth parameter p = 1. 
All values of gm are given in units of a\\huj\\, where a\\ = is the harmonic oscillator length in longitudinal direction. In 

the tilted potential used for the tunneling measurements [Fig. 1(c) of the main text], we calculate oj|| by expanding the combined 
potential K p ti ca i(z) + K iagn etic(z) (Eq. SI and Eq. S2) around its minimum. 


V. Measurement of occupation-number distributions 

To measure the mean occupation numbers on single-particle trap levels, we perform a series of measurements where we spill 
all population above a certain trap level 1 from the trap (Section II) and count the number of remaining atoms (Section III). The 
mean number of remaining atoms corresponds to the sum of the populations on trap levels 0 to i. For the ground state, the 
mean occupancy is directly given by N^°\ For excited states, we obtain the mean occupancies by subtracting from iVW. 















8 



0 200 400 600 800 1000 1200 1400 

Magnetic Field [G] 


FIG. S2: Interaction strength of a quasi-one-dimensional system. ID coupling constant gm between atoms of the two lowest hyperfine states 
of 6 Li as a function of the offset magnetic field. At a magnetic field of approximately 779 G, the 3D scattering length a 3 D equals \/2a±/C 
and g id diverges in a confinement-induced resonance (Eq. S4) 


(a) (b) 



FIG. S3: Correction of occupation-number distributions. Corrected (blue circles) and uncorrected (black squares) occupation probabilities of 
the spin-down atom on single-particle trap levels for a (2,1) system (a) and a (3,1) system (b). For details on the measurement see Fig. 4 of 
the main text. Error bars denote the la statistical uncertainties. 


When measuring the occupation-number distribution of a single spin-down atom in the trap (Fig. 4 of the main text), we 
correct for the finite fidelity of our experiment. During each experimental run, we measure the probability P zero to detect zero 
instead of one spin-down atom in a trap that still contains many trap levels. In an ideal experiment, we would expect this 
probability to be zero. Due to the finite fidelity of both the preparation of the initial state and the detection of the atom number, 
we measure values of P zer0 ss 10 % in all experimental runs. We divide each measured mean atom number TVM by (1 — P zel0 ) 
to correct the data for this finite experimental fidelity. Figure S3 shows the effect of this correction for the data of Fig. 4 of the 
main text. 


VI. Confirmation that final state is ferromagnetic 

Our tunneling model predicts that above the CIR tunneling into the FM final states becomes predominant. To confirm this 
prediction, we show that both the initial (2,1) and (3,1) systems end up in a FM (TV — l)-particle state after the tunneling 
of one particle deep in the super-Tonks regime. For the initial (3,1) [(2,1)] system, the tunneling process is performed at an 
interactions strength of —1/gm = 0.553 ± 0.015 [— 1/giu = 0.495 ± 0.015] and a trap depth parameter p = 0.846 ± 0.025 
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[p = 0.759 ± 0.023]. We ramp the final (N — l)-particle system back through the fermionization regime to zero interaction 
strength. Then, we measure the occupation number distributions on the few lowest single-particle trap levels (Section V) as 
shown in Fig. S4. Within our experimental fidelity, both distributions coincide with the expectation for a FM state. Note that in 
this measurement, the trap for the (3,1) system was significantly deeper than in the measurements for Fig. 3(b) of the main text 
(See Section VII). This could explain the deviation between the data of Fig. 3(b) of the main text and the expected limiting value 
of = 0.75 in the super-Tonks regime. 

(a) (b) 




Mean Occupancy Mean Occupancy 

FIG. S4: Ferromagnetic final states. Combined occupation-number distribution of both spin-up and spin-down atoms on single-particle trap 
levels after tunneling of one atom in the super-Tonks regime, (a), The final two-atom state of an initial (2,1) system after tunneling of one atom 
at — 1/<?id = 0.495 ±0.015. (b), The final three-atom state of an initial (3,1) system after tunneling of one atom at — 1/gw = 0.553 ± 0.015. 
In both cases, the occupation-number distribution was measured after ramping the coupling constant back to gw = 0. The red and gray 
symbols show the occupation-number distributions of the noninteracting states that are adiabatically connected to eigenstates of the ground- 
state multiplet in the spin-chain regime. In particular the red symbols correspond to the expectation for a ferromagnetic state. Error bars 
denoting the lcr statistical uncertainties are smaller than the symbols. 


The spatial wave function of a system of fermions in a ferromagnetic state is completely antisymmetric, which prohibits any 
two particles to be at the same place. Hence, a system of ultracold fermionic atoms with contact interactions in a FM state has to 
be noninteracting. We confirm this for the final two-particle state after letting one atom tunnel from the (2,1) system above the 
CIR by measuring that independent of the interaction strength always one atom occupies the ground state of the trap (Fig. S5). 



FIG. S5: Confirmation that the FM state is noninteracting. Ground-state population of the final state after tunneling of one atom from a (2,1) 
system at — 1/gw = 0.495 ± 0.015 as a function of magnetic field. The labels show the inverse ID coupling constant —1 /gw at each point in 
units of [a||ftu;||] -1 . The small increase of the ground-state population at 600 G and 400 G can be explained by the reduced magnetic moment 
of 6 Li at low magnetic fields, leading to a finite probability of detecting the population of the first excited state. Error bars denoting the la 
statistical uncertainties are smaller than the symbols. 
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(b) 
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FIG. S6: Trap parameters during tunneling measurement. Magnetic fields, initial and final trap depth parameters, ramp times, spill times, and 
inverse ID coupling constants for the measurement of spin-down tunneling (Fig. 3 of main text) for an initial (2,1) system (a), a (3,1) system 
(b), or a (2, 2) system (c). The final trap depth parameter pfi na i during the tunneling process is optimized to let exactly one atom [in (a), and 
(b)] or two atoms [in (c)] tunnel out of the trap. Each row corresponds to one data point in Fig. 3 of the main text. 


VII. Measurement of probabilities of spin-down tunneling 

To measure the probability of spin-down tunneling (Fig. 3 of the main text), we let a specific number of atoms (one atom 
in the (2,1) and the (3,1) system, two atoms in the (2, 2) system) tunnel from the trap and measure the probability of spin- 
down tunneling. To do this, we superimpose the optical dipole trap (Eq. SI) with a magnetic field gradient (Eq. S2) to tilt the 
trap. Thereby we create a potential barrier between the in-trap states and the continuum of states outside the trap. To initiate 
the tunneling process, we ramp down the intensity of the trapping light within a time f ramp and thereby lower the trap depth 
parameter p (Eq. SI) from its initial value pj n j t iai to a final value pg na i- After a time f sp jn, during which atoms can tunnel out of 
the trap, we ramp the power back to its original value to stop the tunneling process. In Table S6a, b, and c the parameters for 
the tunneling measurements on the (2,1), the (3,1), and the (2, 2) system are listed. For the (2,1) and the (3,1) systems, the 
parameters lead to the tunneling of exactly one atom, while for the (2,2) system two atoms tunnel out of the trap. 

Throughout the experimental runs, we perform control measurements in which we detect the number of atoms ( N') directly 
after the tunneling process (see Fig. S7(a), (b), and (c) for the control measurements on the (2,1), the (3,1), and the (2, 2) 
system). The blue circles show the probability to measure the expected N' (2 atoms for the (2,1) system, 3 atoms for the (3,1) 
system, and 2 atoms for the (2,2) system). This probability is always smaller than one due to the finite experimental fidelity of 
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FIG. S7: Raw data of spin-down tunneling. Probabilities to measure N' atoms [(a), (b), (c)] or Ni spin-up atoms [(d), (e), (f)] after the 
tunneling of one atom [(a), (b), (d), (e)], or of two atoms [(c), (f)] at different interaction strenghts. (a) and (d) are for an initial (2,1) system, 
(b) and (e) are for an initial (3,1) system, (c) and (f) are for an initial (2, 2) system. Error bars denote the ler statistical uncertainties. 


our measurements and due to three and four-body losses. Since our experimental fidelity is on the order of 90%, we can assume 
that shots with two or more atoms missing [black squares in Fig. S7(a), (b), (c)] are primarily a result of three and four-body 
losses. For the (2,1) system, three-body losses mostly appear for tunneling close to the resonance [black squares in Fig. S7(a)], 
while for the (3,1) and the (2, 2) systems trap losses can be observed for tunneling at all interaction strengths [black squares in 
Fig. S7(b), (c)]. 

After the initial tunneling process, we selectively remove the spin-down component from the trap and measure the number of 
spin-up atoms (NL) in the trap [Fig. S7(d), (e), and (f) for the (2,1), the (3,1) and the (2,2) systems]. In these figures, the blue 
circles correspond to the probability of spin-down tunneling. 

We correct the probability of spin-down tunneling for three and four-body losses, by subtracting the probabilities of trap losses 
[black squares in Fig. S7(a), (b), (c)] from the respective probabilities to detect 0 spin-up atoms (0, or 1 spin-up atoms in the 
(3,1) system) [black squares in Fig. S7(d), (e), (f)]. Since the sum of all probabilities has to be normalized, this correction 
increases the values of both the red and blue symbols in Fig. S7(d), (e), and (f), while keeping their respective ratio constant. 
The corrected probabilities are shown in Fig. S8. The blue circles in this figure are the corrected probabilities of spin-down 
tunneling that are also shown in Fig. 3 of the main text. 


VIII. Theoretical model 
Spin-chain model 

In the fermionization regime of diverging interaction strength, a ID multicomponent system of N fermionic atoms acquires 
the density distribution of N identical noninteracting fermions. In this limit, the ordering of the atoms is fixed and a one-to-one 
correspondence between the full wave function of space and spin and a pure spin wave function can be established [12, 20]. This 
implies that all observables can be mapped from the full continuous Hilbert space into a discrete spin space and vice versa. This 
mapping therefore allows to treat a ID system in the fermionization regime with a spin-chain model and therefore constitutes a 
substantial simplification of its theoretical description. 
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FIG. S8: Corrected data of spin-down tunneling. Corrected probabilities of the number of spin-up atoms Nl after the tunneling of one atom 
[(a), (b)] or of two atoms (c) at different interaction strengths, (a) is for an initial (2,1) system, (b) is for an initial (3,1) system, (c) is for an 
initial (2, 2) system. The blue circles in (a) and (b) correspond to the probability of tunneling of the spin-down atom. The blue circles in (c) 
correspond to the probability of tunneling of both spin-down atoms. Error bars denote the lcr statistical uncertainties. 


For a spin state |x) consisting of N spin-1/2 particles, we will use the notation 


lx) = |mi) <g>... <g> \m N ) = |mi,..., m N ) , (S5) 

where the spin state | rrii) of the ith particle in the spin chain can either be up (|f)) or down (| ],)). 

The spin-chain Hamiltonian of N spins is given by [12] 


H (N) 



N— 1 \ N -1 

Ji J 1 + JiPi,i+l- 


i=l 


i—1 


(S6) 


Here, is the energy of N spinless noninteracting fermions in ID, i refers to the zth particle, P tJ + i permutes neighboring 
spins and is the exchange coupling between neighboring spins. Ji is given by the formula 


N\h A 

m 2 gm 


dzi■ ■ ■ dz]y5(zi 


z i+ i)0(zi, ...,z N ) 


dipF 

dzi 


(SI) 


Here, tpp is the ground state of N spinless noninteracting fermions in ID and the function 0(z\ ..... z.v) is 1 if zi < • • • < zn 
and zero otherwise. It is important to note that Ji is proportional to the inverse interaction strength (1/<7 id) and that it is 
approximately proportional to the local density cubed and therefore depends on the trap geometry. 

To determine the exchange couplings Ji for our tunneling measurements, we calculate the quasi-bound states ipp in the tilted 
trap (Fig. lc of main text) as follows: First, we approximate the potential by the Taylor series around its local minimum. We 
adjust the maximum order of the Taylor series such that the potential is well approximated up to the barrier maximum and that 
the Taylor series is larger than the barrier maximum outside the quasi-bound region. Then, we calculate the spectrum and the 
eigenfunctions of the approximate potential using the eigenfunctions of the harmonic oscillator at the local minimum. 


Tunneling model 

In our model, we consider tunneling of one atom from the trap. Hence, the model connects an initial spin state |i), which is 
an eigenstate of Hi N \ with a final spin state \ f, t) = |/) ( 8 ) |t). In the final state, |/) is an eigenstate of Hi N 11 and describes 
the N — 1 particles that remain in the trap and |t) indicates the spin orientation of the particle that left the trap. 

The energy E, j of the tunneling particle is equal to the difference between the energies of the initial and final in-trap states. 
Close to the CIR, the energy shifts of all states are linear in — 1/<7id (Fig. 2 of the main text) and the eigenenergies of \i) and |/) 
are well approximated by the solutions of the spin-chain model. In this limit, their difference is 

Eij = (i\H^\i) - (S8) 

Further away from the CIR, we interpolate between the linear energy shifts of the spin-chain model and the corresponding 
saturation energies at —1/gw = ±oo using the rescaled energy shifts of the (1,1) system in the ground state [38]. 
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The weight of each tunneling channel Tjj for tunneling from |«) to |/, t) is approximated by 

Tij oc \(i\f,t)\ 2 E i je~ 2 ' y ^ EiJ \ (S9) 

That means that it is proportional to the squared spin overlap (|(i|/,f)| 2 ), the energy of the tunneling particle (E, : f) and the 
probability (e -27 ^-^) that a particle with energy Eij tunnels through the barrier. The energy-dependent tunneling parameter 
7 (Eij) is given by 

7 (Eij) = \jj dz\J2m [V (z) — E i:f ] . (S10) 

Here, z\ and z 2 are the two largest solutions of V(z) — Eij = 0. 

The probability to tunnel from state |i) to state |/) is given by 


Pi,f = 


Ti 


i,f 




(Sll) 


where the sum is over all possible final states \f). 

Exactly on resonance, — 1/ <j \□ = 0, the energy differences Eij are equal for all final states (Fig. 2 of the main text). Therefore, 
the energy-dependent term Eije~ 2 ' y( ' EiJ 1 in Eq. S9 is equal for all tunneling channels and can be canceled from Eq. Sll. In this 
case the probability to tunnel from state | i) to state \f,t) simplifies to Pij = \(i\f,t}\ 2 since |(*|/ , ,i , )| 2 = 1- The limiting 
behavior of P t j far left and right from the CIR is determined by the predominant tunneling channels of these regimes. For small 
repulsive interactions (left blue arrow in Fig. 2 of the main text), tunneling into the antiferromagnetic (N — l)-particle in-trap 
state is predominant since its energy is much smaller than the energies of all other (N — l)-particle states. For small attractive 
interactions (right blue arrow in Fig. 2 of the main text), tunneling into the FM (N — l)-particle in-trap state is predominant 
since the ordering of the energy levels is inverted. This is also confirmed by the measurements of Section VI. 


IX. Calculation of spin-down tunneling probabilities 

Using Eq. Sll, we calculate the probabilities for spin-down tunneling P±(—l/g\ D ) as a function of interaction strength. Here, 
P^(—1/<7 id) is the probability of tunneling from an initial state |z) into the final state |/) (g) 1 1), where the in-trap state |/) 
only contains spin-up atoms and t) only contains spin-down atoms. The initial and final in-trap states (|z) and |/)) of our 
tunneling model are always eigenstates of the spin-chain Hamiltonian (Eq. S6). In particular, the initial state will always be the 
antiferromagnetic ground state | AFM) of the repulsive spin-chain model with N atoms. 

We will solve the (2,1) and (2,2) systems analytically to make the physics more transparent. To achieve this, we will first 
determine the eigenstates of the (2,1) system in a symmetric trap and then show how they mix in a tilted trap (See Fig. 1(c) of 
the main text and Section VII). In the (2,2) case, we will first determine the two eigenstates of S 2 with 5 = 0. Again, we will 
calculate the spin-chain Hamiltonian within the subspace of these two states and show how they couple in a tilted trap. We will 
show that both systems are closely related, since the (2, 2) system ends up in the AFM states of the three-particle system with 
almost 100% probability after the first tunneling process. 


(2,1) system 

Within the spin basis |1) = | f, |2) = | tst)» anc l |3) = | 4-, tit) the spin-chain Hamiltonian reads 

(-■h J2 0 \ 
i/f) = Ef II + I J 2 -J x -J 2 Ji ■ 

V 0 Jr -Jj 

It has in a symmetric trap (Jr = J 2 ) the eigenstates [12] 

I AFM 3 ) = -J=(| t,t,t> - 2| tit,t) + I t,t,t)), 

|/M 3 ) = -L(|t,t,t)-It,t,t)), 


(S12) 


(S13) 


(S14) 
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I fm 3 ) = (| t, t, 4.) +11,4, t) + 14, t, t)) 


(S15) 


(3) (3) (3) (3) (3) (3) 

and the eigenenergies E AFM = E F — 3Ji, = E F — J\, E FM = E F . In this eigenbasis of the symmetric trap, the 
spin-chain Hamiltonian reads 


H < 3) = E { F ] 1 + 


—1( Ji + J 2 ) 
#(Jr- J 2 ) 


0 


#(Jr-J 2 ) 0 
— \{J\ + J 2 ) 0 
0 0 


fS 16) 


Therefore, an imbalance of J\ and J 2 coherently mixes the AFM with the IM state, while the FM state is decoupled. This can 
be understood, since both the AFM and the IM state have an absolute value of the total spin of S = 1 /2, while the FM state has 
S = 3/2. The eigenstates of in a nonharmonic trap can hence be written as 

| AFM') = cos 11 AFM 3 ) - sin 11 IM 3 ), (S17) 


rv ry 

| IM' 3 ) = cos -\IM 3 ) + sin -I AFM 3 ), (S18) 

and its eigenenergies are E { A l M , = E { F ] - (Ji + J 2 ) (l + 3 ^;), E lm' = E f ] ~ ( J i + h) l 1 “ 2 ^) with the mixin g 
angle a = arctan j 4 ^)- 

The (2,1) system is initially prepared in state \AFM 3 ). Exactly on resonance, —1/gw = 0, the ordering of the particles is 
fixed by their mutual repulsion such that only the rightmost particle can leave the trap (Fig. 1(c) of the main text). The probability 
for spin-down tunneling, Pj_, is hence just the probability of the rightmost spin of the \AFM' 3 ) state to point downwards, 

n(-l/ 5 1D = 0) = |(t,t4|^M')| 2 = (cos|-^-sin|^ . (S19) 

In the symmetric trap (a = 0), the probability for spin-down tunneling is P^{—l/gw = 0) = 1/6 = 16.7%. In the tilted trap, 
q oc Ji — J 2 is positive, since the density between the left particle and the particle in the middle is larger than the density between 
the particle in the middle and the right particle (Fig. 1(c) of the main text). P^{—l/gw = 0) is hence lower than 16.7% in the 
tilted trap. We obtain for our experimental trap parameters during the tunneling process a mixing angle of a ~ 20° resulting in 
P\.{—l/gm = 0) ~ 8 % for the | AFM/) state (Fig. 3(a) of the main text). 

Away from resonance, — 1 /(?id 7 ^ 0, one has to calculate the weights of all the possible tunneling channels. The eigenstates 

of #i 2) and ^ 3) are |FM 2)1 ) = | t,t4), I™ 2 , 2 > = ^(| t4) + | 4,t»l t>. and |AFM 2 ) = ^(| f,l) - I l,t»| t>- These 

( 2 ) ( 2 ) ( 2 ) ( 2 ) 

eigenstates have the eigenenergies E FM = E F and E AFM = E F — 2J\ (now calculated with J\ of the two-particle in-trap 
system). The probability for spin-down tunneling, P±, is the probability to tunnel from the initial state \ AFM//) into the final 
state |FM 2 ,i) = | t> 4-)- Far left from the CIR, tunneling into \AFM 2 ) is predominant and one obtains the limiting value 

P^-l/gw -> - 00 ) = 0 . (S 20 ) 


Far right from the CIR, tunneling into \FM 21 ) and FA/ 2 . 2 ) is predominant and one obtains the limiting value 


P (_M t = _ \(FM 2t i\AFM 3 )\ 2 _ 2 

U /9m ' \(FM 2A \AFM 3 )\ 2 + \(FM 2 , 2 \AFM 3 )\ 2 3 


(S21) 


(see Fig. 3(a) of the main text). Interestingly, this limiting value is independent of the mixing angle a and hence the tilting of 
the trap. In between, P± = P\afm^),\-\-,^) is calculated using the parameters pfj na 1 and — 1/gw of Table S6(a). This leads to the 
red solid curve in Fig. 3(a) of the main text. Similar calculations for the IM and FM initial states lead to the dotted and dashed 
curves of Fig. 3(a) of the main text, respectively. 

We assumed for the above calculations that the trap is tilted so slowly that the spin configuration can follow the change of J\ 
and ,/ 2 induced by the change of the total density. This assumption is wrong close to the CIR, where the spin configuration of the 
initially deeper trap may be frozen. We would hence expect a larger P / close to the CIR. A full dynamical calculation showed, 
however, that this effect is not sufficiently strong to explain the resonance feature of P / close to the CIR. 
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(3,1) system 


Within the spin basis 11) = 
reads 


12) = | 13) = | and |4) = | the spin-chain Hamiltonian 


HW = E { ph + 


(~J 3 j 3 0 0 \ 

■h —Ji — J 3 -h 0 
0 J 2 — Ji — J2 Ji 
V 0 0 Jr ~.h) 


(S22) 


The (3,1) system is initially prepared in \ AFM 4 ), which is the ground state of this Hamiltonian for g > 0. We determine 
this state numerically. We have to consider tunneling into the final states |FM 3j1 ) = | |PM 3j2 ) = | FM 3 }\ f), 

|/Mg') = \IM 3 )\ t), and \AFM 3 ) = |AFM 3 )| f), where \FM 3 ) is given by Eq. (S15), \IM 3 ) is given by Eq. (S18), and 
\AFM 3 ) is given by Eq. (S17). Exactly on resonance, the probability for spin-down tunneling is just given by 

n(-i/siD = 0) = \AFM a )\ 2 . (S23) 


We obtain for the harmonic trap a value of P^(— 1/giD = 0) = 5.1% [12], This value is again lowered by the tilting and we 
obtain for the experimental parameters at the CIR (see Table S6(b)) a value of P^(— 1/gw = 0) ~ 1%. Far left from the CIR, 
tunneling into | AFM'^) is predominant and one obtains again the limiting value P^(—1/giD —> —00) = 0. Far right from the 
CIR, tunneling into the two FM final states, |PM 3 j i) and \FM 3y2 ), is predominant and one obtains the limiting value 


p,_ Un . ,» \(FM 3 A \AFM A )\ 2 3 

U 9m \(FM 3 } 1 \AFM a )\ 2 + \(FM 3 i 2 \AFM 4 )\ 2 4’ 


(S24) 


which is again independent of the tilting. The experimentally measured value for Pj, is, however, only « 50% in the super- 
Tonks regime (see Fig. 3(b) of the main text). The deviation from the theoretical prediction (« 70%) is presumably caused by 
a stronger contribution of the IM tunneling channel in the experiment. Indeed, the IM tunneling channel contributes stronger if 
the lowering of the potential barrier is assumed to be nonadiabatic, as outlined in the following paragraphs. 

In the fermionization regime, the energy of the energetically highest particle is given by the A 7 th trap level. This energy grows 
linearly with — 1/gio in the fermionization regime and it saturates in the weakly-interacting regimes, — 1 /g m = ±00, at values, 
which are again given by the (quasi-bound) trap levels [24]. The saturated energy shifts can hence never exceed the barrier 
maximum. This situation applies to an adiabatic lowering of the potential barrier before the tunneling process. This condition 
is fulfilled, if the tilted trap is deep enough or if the barrier is lowered slowly when the corresponding trap level approaches the 
barrier maximum. 

If the barrier height is nonadiabatically lowered, the saturated energy shifts may exceed the barrier maximum. In this case the 
saturated shifts are multiples of the level spacing at the Fermi edge and they are not limited by the barrier maximum. The energy 
of the tunneling particle may now exceed the barrier maximum for the energetically most favorable tunneling channels. In that 
case the particle tunnels out with a probability of 1. Moreover, the tunneling energies of the intermediate channels may also be 
slightly below or above the barrier maximum. This is the case for the IM tunneling channel of the (3,1) system in the super-Tonks 
regime. As a result, the weight of the IM tunneling channel is much stronger than in the adiabatic approximation, which leads to 
a substantial reduction of Pj, in the super-Tonks regime. We obtain P^(— l/gm = 0.35) = 38% and P^(— 1/<?id = 0.5) = 61% 
for the (3,1) system in the super-Tonks regime, which is in much better agreement with the experimental data. We finally note 
that, in contrast to the (3,1) system, both models lead to nearly the same results for the (2,1) and (2,2) systems. 


(2,2) system 


Within the spin basis |1) = | t.t.I.I). 1 2) = | |3) = | 14) = | 1,1, t,t). ]5) = 

16) = | the spin-chain Hamiltonian reads 


and 


P) 4 ) =E { p ) 1 + 


(-J2 J-2 0 0 

J 2 —Ji — J 2 — J 3 J 1 0 

0 J\ — J\ — J 3 0 

0 0 0 -J 2 

0 0 J 3 

V 0 j 3 0 


0 

0 

J 3 

J2 


0 

J 3 
0 
0 

J 2 —Jl — J 2 — J 3 J\ 

0 Jl -Ji-J 3 ) 


(S25) 
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The (2,2) system is initially prepared in \AFM a ), which is the ground state of this Hamiltonian for g > 0. It can be determined 
analytically since only the ground state and another excited have an absolute value of the total spin of S = 0. The two eigenstates 
with S' = 0 are 


and 


|S = 0,l) = i(|l)-|3) + |4)-|6)) 

15 = 0) 2 > = 7TI ( |1} " 2|2> +13> +14) _ 2|5) + |6) )' 

In the subspace of these two states, the spin-chain Hamiltonian attains the form 


H^ 


s =o 


= E ( p ] 1 + 


— |(Ji + Ji + J 3 ) — J2 + J 3 ) 

^r{Ji — J2 + S3) —|(Ji + J2 + S3) 


The AFM ground state is hence 

| AFM') = -sin^|S = 0,l)+cos||S = 0,2) 


•s/12 


with the mixing angle /3 = arctan 


cos^ - 4 sin - | | t,t44) - 4s C0S sl t,7t4) + 


73 

3 anc * *^ le enel 'gy E^ 


7rI C0S ? + 2 sin f 


(526) 

(527) 

(528) 

(529) 

(530) 


AFM 1 


= E 


(4) 


resonance, the probability for the tunneling of two spin-down particles is therefore given by 


— (Si + S 2 + S 3 ) ^1 + 2 c * s p ^ • Exactly on 


Pd-i/giD = 0 ) = |(t,t44 \afm’ a )\ 2 = ( 


1 P 1 . P 

—— cos — — - sin — 

712 2 2 2 


(S31) 


We obtain for the experimental parameters at the CIR (see Table S6(c)) an angle of j3 « 21° and hence P±(— I/^id = 0) ~ 4%. 
Away from resonance, we have to consider the sequential tunneling of two particles from the initial \AFM' 4 ) state. In the 
following, we discuss only the case that a spin-down particle has tunnelled in the first tunneling process (the spin-up case is 
analogous). Therefore, we have to consider the tunneling into the final states FjM() = |FM 3 )| 4 ,), = \IM^)\ 4 ,), and 

\AFM” 1 ) = | AFM') | 4 .), where IFM 3 ) is given by Eq. (S15), |FMg) is given by Eq. (SI 8), and |AFMg) is given by Eq. (S17). 
It follows from Eqs. (S15) and (S30) that the spin overlap of \AFM a ) with FMg) is exactly zero independent of the tilting of 
the trap. Moreover, the squared spin overlap with the IM state. 


\(IM^'\AFM a )\ 2 = ^ sin 2 


reaches at most 0.1%, while the squared spin overlap with the AFM state, 

1 


|<AFM'"|AFM7/ = -cos 2 


a — /? 


(S32) 


(S33) 


is at least 49.9% (the remaining 50% belong to the spin-up case). Therefore, even in the super-Tonks regime, where tunneling 
into the IM state is energetically favored, the particle tunnels with at most 1% into the IM state due to the small spin overlap. 
We conclude that the system ends up in the AFM three-particle in-trap state with sa 50% after the first tunneling process. The 
second tunneling process is then the tunneling from the AFM state of a (2,1) system, which has been discussed above. It follows 
that the probability for spin-down tunneling, Pg = P\afm' a ), * s f° r the (2,2) system half as large as for the (2,1) system in 
the whole interaction regime. This is indeed the case, as can be seen in Fig. 3(c) of the main text. 


X. Exact diagonalization 


The occupation-number distributions of Fig. 4 of the main text have been calculated by means of an exact diagonalization of 
the full continuous model. The second-quantized Hamiltonian of the full model reads 


H = 


Z^ 

a=4„f 


j dz¥ ^ {-L^ + +5id . 


dzdz' 4>^(z)4>^(z')8{z — z')'i'.|.(~ , )'^T( 2: ) (S34) 
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with the fermionic field operators 'f' 0 .(z) and the angular frequency oj| of the harmonic oscillator at the local minimum. We 
expand the field operators in the energy eigenbasis of the harmonic oscillator, \E , cr (z) = ]T\ 4>i{z)di t<J . Here, i refers to the ?'th 
energy level, <j>i(z) is the /th energy eigenstate of the harmonic oscillator, and a,. n annihilates a fermion in energy eigenstate i 
and spin state er. The Hamiltonian is then numerically calculated and diagonalized in the many-body energy eigenbasis of the 
noninteracting problem. Here, the finite size of the computer memory restricts the calculation to all noninteracting eigenstates 
below a certain cutoff energy. This allows for an accurate numerical calculation of the low-energy eigenstates of the interacting 
few-body system (for not too many particles and not too strong interactions, |c/ 1D | < 20). The mean occupancy of single-particle 
level (i, a) is the expectation value (aj a d i jCr ) of the corresponding eigenstate of the interacting system. 



